THE EXPONENTIAL MAP AT A CUSP SINGULARITY 



VINCENT GRANDJEAN AND DANIEL GRIESER 

Abstract. We study spaces with a cusp singularity embedded in a smooth 
Riemannian manifold and analyze the geodesies in these spaces which start at 
the singularity. This provides a basis for understanding the intrinsic geometry 
of such spaces near the singularity. We show that these geodesies combine 
to naturally define an exponential map based at the singularity, but that the 
behavior of this map can deviate strongly from the behavior of the exponential 
map based at a smooth point or at a conic singularity: While it is always 
surjective near the singularity, it may be discontinuous and non-injective on 
any neighborhood of the singularity. The precise behavior of the exponential 
map is determined by a function on the link of the singularity which is an 
invariant - essentially the only boundary invariant - of the induced metric. 
The results are proved in the more general natural setting of manifolds with 
boundary carrying a so-called cusp metric. 
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1. Introduction 

Geodesies are among the fundamental objects of differential geometry. On a 
smooth Riemannian manifold the geodesies starting at a point p are classified by 
and smoothly depend on their initial velocity vector, and combine to define the 
exponential map based at p, which in turn yields normal coordinates and important 
geometric information about the manifold. Also, geodesies arise in the study of the 
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propagation of waves on the manifold as paths along which singularities of solutions 
of the linear wave equation travel. 

Much less is known about the behavior of geodesies on singular spaces. Here by 
a geodesic we always mean a locally shortest curve. Our general aim is to analyze 
the full local asymptotic behavior of the family of geodesies reaching, or starting at, 
a singularity. Previously, Bernig and Lytchak BL obtained first order information 
for geodesies on general real algebraic sets X C R™, by showing that any geodesic 
reaching a singular point p of X in finite time must have a limit direction at p. In 
the special situation of isolated conic singularities p, Melrose and Wunsch [MW 
obtain full asymptotic information by showing that the geodesies starting at p define 
a smooth foliation of a neighborhood of p and thus may be combined to define a 
smooth exponential map based at p, cf. Theorem 1 1 . 2| below. 

In this paper we analyze the family of geodesies starting at an isolated cusp sin- 
gularity, defined below. Our method is based on the geodesic differential equations 
and yields full asymptotic information about the geodesies. We will see that there is 
a much richer range of possible local behavior than in the case of conic singularities; 
for example, there is a natural notion of exponential map based at the singularity, 
but this map may be non-injective or discontinuous on any neighborhood of the 
singularity. 

A natural setting for our investigation is the notion of cusp manifold. This 
is a smooth (that is, C°°) manifold, M, with compact boundary, equipped with 
a semi-Riemannian metric, g, which is Riemannian in the interior M and in a 
neighborhood U of the boundary dM can be written 

(1) g\ v - {l-2r 2 S + 0{r z ))dr 2 + r i h 

where r is a boundary defining function for M (that is, r £ C°°{M) is positive in 
M, vanishes on dM and satisfies dr p ^ for all p £ dM), S is a smooth function 
on dM and h is a smooth symmetric two-tensor on U whose restriction to dM, 
denoted go a/, is positive definite. We call g a cusp metric. More precisely, this is 
a cusp metric of order two. Cusp metrics of higher order are defined in Section [8j 
where we also show how our results carry over to this more general case. A given 
cusp metric defines invariantly the quantities 

S £ C°°(dM) and go a/, a Riemannian metric on dM 

(more precisely, S is only determined up to an additive constant) and fixes r to 
third order. For simplicity we will fix a boundary defining function r throughout. 

Cusp manifolds arise as resolutions of spaces X with isolated cusp singularity 
p £ X (of order two) . These are subsets X C K™ so that X \ p is a submanifold of 
R" and so that X is given by the following local model in a neighborhood U' of p: 
Let p = 0. Then, with K + = [0, oo), 

(2) XnU'=P(X), (3 : xl + -> E™, (u, z) H> (z 2 u, z) 

where X C W 1 ^ 1 x R + is a smooth manifold with boundary, whose boundary is 
contained in E = R™ -1 x {0} and so that X intersects E transversally. We call 
X the resolution of X. See Figure [l] Now if g is any smooth Riemannian metric 
on the ambient space W 1 then the restriction of g to the manifold X \ p pulls back 
under /3 to a Riemannian metric on the interior of X, which extends smoothly to 
the boundary. The extension is a cusp metric on X. In this sense the resolution of a 
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space with cusp singularity is a cusp manifold. See Section[2]for a discussion of cusp 
manifolds, cusp singularities and the meaning of S and ggM for cusp singularities. 

Let (M, g) be a cusp manifold. Since g is Riemannian over the interior M, 
locally shortest curves in M are precisely the unit speed solutions of the geodesic 
differential equations. We call these solutions geodesies on M. We are interested 
in those geodesies starting at the boundary, that is, geodesies 7 : (T_,T + ) — > M for 
which lim r _j.T_ K7( T )) = 0- Here T_ glU {—00}. If lim r _yT_ 7(7") exists then we 
call it the starting point of 7. 

In our first theorem we show that, under certain assumptions, geodesies starting 
at the boundary have a starting point on dM, and characterize the possible starting 
points. Recall that a critical point of S G C°°(dM) is a point p G dM where the 
differential dS\ p vanishes. 

Theorem 1.1. Let (M, g) be a cusp manifold. Assume that either S is constant 
or the critical points of S are isolated. 

Then there is r > so that the following is true: Let 7 : (T_,T + ) — » M be a 
maximal geodesic and assume that for some T £ (TL,T + ) the initial part J\(t_,t) 
is contained in {r < ro}. Then T_ is finite and the starting point lim T _s.T_ 7(t) 
exists. It is a point on dM which is a critical point of S. 

Conversely, each critical point of S is the starting point of a geodesic on M. 

Our next theorem says that for constant S we have a smooth foliation by 
geodesies of a neighborhood of p, analogous to the Theorem of Melrose and Wunsch 
[MWj in the conic case. 

Theorem 1.2. Let (M, g) be a cusp manifold. Assume S is constant. Then to each 
p G dM there is a unique geodesic 7 P starting at p at time r = 0. Furthermore, 
there is To > and a neighborhood U of dM in M such that the exponential map 

exp aM : dM x [0, r ) -> U C M, (p, r) h-> 7p (r) 

is a diffeomorphism. 

See the left picture in Figure [2j For a general cusp manifold, where S need not 
be constant, it seems clear that a full neighborhood of the boundary is covered by 
geodesies starting at the boundary, since for any interior point of M any shortest 
curve from this point to the boundary must be a geodesic. We will not attempt 
to make this argument precise here (that is, prove existence of a minimizcr) but 
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Figure 2. Geodesies in the cases S = const, S w < 9/4 and 
<SW(Pmin) > 9/4. The bold line is dM. The points p max and p m i n 
are a maximum resp. minimum of S. 



rather analyze the case where S is a Morse function much more precisely, mostly 
in the case of surfaces. 

For our next theorems we assume that S is a Morse function. That is, all critical 
points are non-degenerate (i.e. the Hessian of S at these points is a non-degenerate 
quadratic form). This implies that the critical points are isolated. 

The behavior of geodesies hitting a non-degenerate critical point p of S can be 



analyzed very precisely in a neighborhood of p. See for example Proposition 4.2 for 
the case of local maxima and minima of S. We will now focus on surfaces M, for 
which we can analyze how these local pictures fit together in a full neighborhood of 
dM. First, we show that there is an exponential map. To simplify the exposition 
we assume in the following theorems that dM is connected. 

Theorem 1.3. Let (Af, g) be a cusp surface with connected boundary. Assume that 
S is a Morse function. There is a parametrization q > jq, q € S 1 , of the geodesies 
starting at dM at time r — so that the map 

(3) ex PaA/ tS 1 x (0,t )^U\8M, ( q , T ) ^ lq {r) 

is defined and surjective for some tq > and some neighborhood U of dM. 

The full asymptotic behavior of exp aM as r — S- can be described explicitly in 
terms of S. 

The parametrization q ^ j q is uniquely determined, up to homeomorphisms of 
S 1 , by the requirement that it preserves cyclic ordering in a suitable sense. However, 
exp aM need not be a diffeomorphism, it may even be discontinuous. See Section [6] 
for details. 

In the next two theorems we give conditions under which more can be said 
about the exponential map. Here and throughout we denote by S^, S vv> the first 
and second derivatives of S with respect to a coordinate if on dM. 

Theorem 1.4. Let (A/, g) be a cusp surface with connected boundary. Assume that 
S is a Morse function, that < | on dM and that S vv never takes the value 
2 at any local minimum, where ip is an arc length parameter on dM. Then the 
exponential map ^ is a homeomorphism for suitable To, U . 

In particular, there is a neighborhood of dM such that U \ dM is foliated by 
geodesies starting at dM. 

The exponential map ^ extends to the boundary by letting exp dM (q, 0) be 
the starting point of 7 g . However, unlike in the case of constant S, the extension 
S 1 x [0,To) U is not continuous in q; this is clear since the image of S 1 x {0} 
is a finite set with at least two elements, a maximum and minimum of 5*. In 
the setting of Theorem |1.4| the cusp metric defines naturally a compactification 
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M' of the interior M different from M, and the relation between M and M' is 
similar to a birational map, see Remark |6.7[ The exponential map exp aM extends 
to a homeomorphism from S 1 x [0,To) to a neighborhood of the boundary of M', 
and could be proved by our methods to have mildly higher regularity than just 
continuity. 

is optimal in the following sense. 



The value | in Theorem 



1.4 



Theorem 1.5. Let (M, g) be a cusp surface with connected boundary. Assume that 
S is a Morse function and that S vv> > | at some minimum of S, for an arc length 
parameter ip on dM . Then the exponential map ^ is not injective for any tq > 0. 

That is, in any neighborhood U of dM there are points through which at least 
two geodesies starting at the boundary pass. 

Theorems 1 1 . 4| and 1 1 . 5| are illustrated in Figure [2j middle and right. 



We also show that the conditions on S in Theorem (1.4) are satisfied if M arises 
from a cusp singularity X and if dX is contained in the boundary of a strictly 
convex subset of R™ -1 which contains the origin, and is never doubly tangent to a 



sphere centered at the origin. See Proposition 5.5 



Theorems 1.2 1.3 1.4 and 1.5 may be summarized as follows, in the case of 
surfaces: There is a well-defined exponential map. For constant S it is a smooth 
diffeomorphism near the boundary. If S is not too far from constant then it is a 
homeomorphism, though not on the boundary. If S is far from constant then it 
may be not injective and also discontinuous. 

Main ideas, outline of the proofs. Let g* be the metric on the cotangent 
bundle T* M dual to g. Consider the energy function (Hamiltonian) E = Ag* and 
the associated Hamiltonian vector field W on T*M. Then geodesies on M are 
the projections to M of integral curves of W. Unit speed geodesies correspond to 
integral curves on the energy hypersurface {E = 5}- The degeneracy of g at dM 
implies that E and hence W are undefined over dM. To make this explicit, let 
(r, ip), if — (tp 1, . . . , (fim-i) be local coordinates near a boundary point of M and 
denote by £, rj = (771, ... , J7 m _i) the dual coordinates on the fibers of T*M. Since g 
is a positive definite quadratic form in dr and r 2 dip whose coefficients are smooth 
functions of r, <p up to r = 0, the function E is a positive definite quadratic form in 
£, ^2 with coefficients smooth up to r = 0. Therefore, rescaling 

(4) 0=1, E(r,cp,^e)=E(r,cp,t,r 3 0) 

yields a function E which is smooth up to the boundary r = 0, and simple calcula- 
tions using the specific form |[lj) of g show that the associated Hamiltonian vector 
field W (which is W written in coordinates r, tp, £, 9) is \ times a vector field V 
which is smooth up to the boundary r = and also tangent to the boundary. 

Clearly, V and W have the same integral curves up to time reparametrization, so 
we need to analyze the integral curves of the rescaled vector field V. It is essential 
for our analysis that V is sufficiently non-degenerate at its singular points to make 
a precise analysis of its integral curves possible. For example, the singular points 
of V are hyperbolic if S is a Morse function. 

It may seem strange to use the rescaling Q instead of fi = -pj, which is more 
naturally associated to cusp manifolds and is sufficient to make the energy a smooth 
function. But it turns out that the Hamilton vector field, written in coordinates 
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r, p, £, fj,, needs to be multiplied by r 2 instead of r to yield a smooth vector field, 
and that the resulting vector field is highly degenerate near its singular points, so 
a precise analysis via linearization would not be possible. 

The rescaling Q may be given an invariant description as follows. This is the 
natural setting for our results and methods, but is not strictly necessary to under- 
stand most of the paper. Introducing 9 corresponds to replacing the vector bundle 
T* M by a rescaled cotangent vector bundle which we denote by 3 T*M, defined as 
the unique vector bundle over M whose space of sections is the set of smooth one 
forms a on M satisfying a(V) — 0(r 3 ) for all smooth vector fields V tangent to 
third order to the boundary, i.e. satisfying dr(V) — 0(r 3 )j^] In coordinates, this 
space of sections is spanned by dr, r 3 dp\, . . . ,r 3 dp m _\ over C°°(M). Then E is 
a smooth function on 3 T*M and the rescaled geodesic vector field V is a smooth 
vector field on 3 T*M. 

We also need to describe the relation of 3 T*M to T*dM. The rescaled cotangent 
bundle 3 T*M is a manifold with boundary d{ 3 T*M) = 3 Tq M M , and tt*t (where 
7r : 3 T*M — > M is the projection), denoted r for short, is a boundary defining 
function. The coordinate £ is invariantly defined on 3 Tg M M, and for any £o € K 
the affine subbundle £ = £ 01 3 Tg M M may be naturally identified with T*dM. 
This identification is given in coordinates by £0 dr + 9r 3 d(p 1— > 9 dip, where we write 
9dp = YT=i 6 i dL Pi- That is > in coordinates the embedding T*dM -> 3 Tq M M is 
simply given by (p, 9) i-> (0, tp, £oj &) where on the left 9 plays the role of the fiber 
coordinate on T*dM. 

Our aim is to understand the dynamics of V close to the boundary, so the 
dynamics of the restriction of V to the boundary plays an important role in our 
analysis. Unit speed geodesies leaving the boundary have f = 1 there, so this 
restriction may be considered as a vector field on T*dM under the identification 
just described. We will see that the boundary dynamics is given by a damped 
Hamiltonian system with Hamiltonian 

E d = S+\z* dM on T*dM. 

More precisely, in local coordinates it is given by the equations 

(5) <P = E$, 6=-E°-W. 

Here Eg — (^f~, • ■ ■ , d ® E - ) etc. The Hamiltonian E d corresponds to a particle 
moving on the Riemannian manifold [dM, ggjvf) m the potential S. 

The theorems now follow from a precise analysis of V. Singular points p € 
3 Tq M M of V (that is, points where V vanishes) correspond to critical points p £ 
dM of S via projection, and integral curves of V leaving p correspond to geodesies 
starting at p. These integral curves foliate the unstable manifold M£ c 3 T*M of 
p. If S is constant then each p € dM is a critical point, and V is transversally 
hyperbolic with respect to the submanifold of T*M formed by the corresponding 



points p, and the unstable manifold theorem yields Theorem 1.2 For general S, 
the analysis of the exponential map breaks up into three parts: 

a) Understand the flow on Mp for each critical point p. 



^For simplicity, we assume a fixed boundary defining function r is chosen, although to define 
3 T*M it is enough to require the choice of r up to changes of the form r k-> a±r + aiT 2 + 0(r 3 ) 
where a\ > 0, 0,2 are constants. 
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b) Understand the regularity of M u = Up-^pj that is, how the various M£ 
fit together. 

c) Determine whether M" projects diffeomorphically to M in a neighborhood 
of dM, under the projection 3 T*M -» M. 

If S is a Morse function then the singular points of V are isolated, and the linear 
part of V is invertible at each singular point; this allows to do a) and answer c) 
affirmatively (for Mp*) near each p. The main task now is to understand the global 
boundary dynamics, i.e. the global behavior of the unstable manifolds = dM™ 
of V restricted to the boundary. In the case of surfaces this global analysis can 
be done. This allows to do a), b) and c) in a full neighborhood of the boundary. 
Here some technicalities about linearizations enter; the assumption S vv ^ 2 at 
minima in Theorem 11.41 is a non-resonance condition needed to ensure existence 
of a C 1 linearization. However, we believe that the theorem is true without this 



assumption. From this we then deduce Theorems 1 1 . 1 1 |L3 1.4 and 1.5 



Structure of the paper. In Section[2]we introduce cusp manifolds and show how 
they arise from cusp singularities. In Section [3] we calculate the rescaled geodesic 
vector field V and analyze its singular points, and prove Theorem [L2j In Section [4] 
we analyze the local behavior of the flow of V for a Morse function S, and discuss 
linearizations. In Section [5] we analyze the boundary dynamics, which we use in 
Section [6] to prove the remaining theorems, after defining the exponential map. 
In Section [7] we give some examples and in Section [8] we indicate the adjustments 
needed for higher order cusps. Here the same theorems hold, with essentially the 
same proofs, except that | has to be replaced by ^(k-i) ^ or a CUS P °f order k. 

Further related work. The exponential map may be used to study differential 
geometric quantities near a singularity. For example, the description of the as- 
ymptotic behavior of the exponential map may be used to deduce the asymptotic 
behavior of the volume of balls centered at the singularity, as the radius tends 
to zero. In the case where p is a smooth point, the coefficients in this expansion 
are related to the curvature at p. In }Gri2) such an expansion was derived for 
real analytic isolated surface singularities; however, the balls were defined extrin- 
sically, i.e. with distance defined as distance in the ambient space R™. Also, the 
bi-Lipschitz geometry of (singular) algebraic subsets of R n was investigated exten- 
sively in [Mo[ IPa| lB"H IValj [Va2 . However, very little is known about the precise 
intrinsic geometry of these sets, and the present paper provides a first insight into 
this question. 



2. The geometric setting 

In this section we first define the notion of a cusp singularity of order two of 
a subset X of a smooth manifold Z. Since this is a local notion, we may assume 
Z = R™. We give two equivalent characterizations. Then we introduce the notion 
of cusp metric on a manifold with boundary, which captures essential features of 
the metrics induced on a space with a cusp singularity of order 2 from a smooth 
ambient metric. A similar discussion applies to cusp singularities of higher order. 

First definition of a cusp singularity of order 2. Let X be a subset of M. n . 
The tangent cone C(X,p) of X at p <E X is the set of limits of oriented secant 
half-lines through p from points in X: 
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C(X,p) := {rv : r > 0,u G S"" 1 , 3{ Pk ) k eX\{ P }: Pk ^ P and lim fc = is}. 

Assume A is a subset of M. n so that X \ { P } is a submanifold. We say that X has 
a cwsp singularity of order two at p if the tangent cone C (X, p) is a single half- line 
and the singularity can be resolved by a single quadratic blowing-up. That is, after 
choosing coordinates (x±, . . . , x n -i, z) in E™ so that p = and C(A, 0) = M>o<9 z 
the space X is given as in pi). 

Definition of cusp singularity of order 2 in terms of iterated blowings-up. 

Equivalently, cusp singularities of order 2 can be characterized by an iteration of 
standard point blowings-up. We explain this here since it lends itself to natural 
generalizations and since it allows a simple derivation of the standard form for the 
metric below. 

Recall the notion of (oriented) blowing-up of a manifold Z in a point geZ, which 
is a geometric way of introducing polar coordinates around q. This is a manifold 
with boundary, denoted [Z, q] , together with a smooth map f} q : [Z, q] — > Z (the 
blowing-down map) which maps the boundary d[Z, q] to q and is a diffeomorphism 
from [Z, q] \ d[Z, q] to Z\ {q}, and which is, locally near d[Z, q] resp. q, given by the 
following model: Z = M n , q = 0, and the n [Z, q ] = R+ x S™" 1 and P q (r,u) = raj. 
Thus d[Z,q] = {0} x S"^ 1 . See [Melj or |Gril| for a more in-depth discussion, in 
particular of the coordinate invariance of this notion and of its generalization to 
manifolds with corners. Of this we only need the case where Z is itself a manifold 
with boundary and q € dZ. Then [Z, q] is a manifold with corners; the local model 
is Z = M"- 1 xR + ,q = 0, [Z, q] = R+ x S+~\ where S™" 1 = S"" 1 n (M™ _1 x E+) 
is the upper half sphere and j3 q is as before. In this case, j3 q maps {0} x S" _1 to 
and is a diffeomorphism between the complements of these sets. In either case 
(3~ 1 (q) is called the front face of the blowing-up. 

If X C Z then the strict transform of X under the blowing-up of p is the 
closure of the pre-image of X \ {q}, so /3*A := f3 q 1 (X \ {q}). A p-submanifold of 
a manifold with boundary Z is a submanifold X c Z such that dX C dZ and 
X hits dZ transversally. This extends to manifolds with corners Z\ we only need 
the straightforward case where X intersects the boundary only in the interior of a 
boundary hypersurface. 

Let Z be a manifold. We say that X C Z has a conic singularity at q G X if 
X \ {q} is a submanifold and the strict transform of A is a p-submanifold of [Z, q] . 
If Z is a manifold with boundary and q G (5Z then we require, in addition, that dX 
be contained in the interior of the front face of [Z, q] . 

A subset X <Z Z has a cusp singularity of order 2 at p G A if its strict transform 
A' = (3p(X) C 2f' := [Z, p] intersects dZ' in a single point q and has a conic 
singularity there. Thus, the singularity of A can be resolved by first blowing-up p 
and then q. The resolution is the strict transform 

(6) X — (f3 p o (3 q )*X. 

To see the equivalence of the two definitions of cusp singularity of order 2 recall 
the definition of projective coordinates on [R n ,0]: Let (x,z) be coordinates on 
E™ where x = (x%, . . . , x n _i) € M™" 1 . Then projective coordinates (y,z), y — 
{yi, ... , y n -i) € z G M+ can be introduced on the 'upper half U = {(r, lj) G 

R + x S" 1-1 : w„ > 0} of [K™,0]. These are uniquely determined by requiring that, 
in terms of these coordinates, the blowing-down map is (3 (y,z) = (zy,z). (The 
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relation to the r, w-variables is yi — ^- cjj, z = ru n , but this is rarely needed.) 
There are also projective coordinates covering the remainder of [R n ,0], but we do 
not need them here. 

The two definitions of cusp singularity of order 2 of X C Z = M n at p = 
are related as follows. Let (3q : []R n ,0] — > M. n be the blowing-down map. Points 
of the front face (9[R™,0] correspond to directions at 0, so C(X, 0) = {rv : v E 
d(f3QX),r > 0}. Thus, q in the second definition spans the half line C(X, 0). 
Assume C(X, 0) = IR>o<9 z . Then we may use the projective coordinates introduced 
above on Z' = [M. n , 0] near q, so we identify a subset of Z' by M n_1 x K + , and then 
q = 0. Now blow-up q in Z' and denote the blowing-down map f3 q . Use projective 
coordinates again, now on [Z',q] Z> [K n_1 x R + ,0]. Call these (u, z), so that 
j3 q {u,z) — (zu,z). The fact that, in the second definition, the strict transform of 
X' hits the front face of \Z' , q] in its interior means precisely that, in a neighborhood 
of the boundary, it is contained in the domain of definition of these coordinates. 
Then the combined blowing-down map (/3 P o f3 q )(u,z) = (z 2 u,z) is precisely the 
quadratic blowing-down map f3 in ([2]). 

Cusp metrics. 

Definition 2.1. A cusp metric on a manifold with boundary M is a smooth semi- 
Riemannian metric g on M for which near any boundary point there are coordinates 
(r, ipi, . . . , ip m ^i), with the boundary defined by r = 0, in which g takes the form 

(m— 1 m— 1 \ 

2 b i dr d ^ + c i] d Vt dfj 
i=l i,j=l J 

where S = S(ip), bi — bi(r, ip) andcij = Cij(r, ip) are smooth (Mi<i X}7j'=i c ij(0j d( Pi d( Pj 
is positive definite for each ip, so defines a Riemannian metric on the boundary. 

Observe that the form of metric |7|) is unaffected by a change of the (/^-coordinates. 
However, this is not true for the coordinate (boundary defining function) r. More 
precisely, we have: 

Lemma 2.2. A cusp metric on M determines invariantly the function S on dM , 
up to additive constants, as well as a Riemannian metric gdM on dM, given by 
SdM — X)™j=i c ij'(0' "Z 5 ) dipidcpj in coordinates. The metric also fixes the boundary 
defining function r up to replacing it by r + br 3 + 0(r 4 ) 7 where b is constant on 
each connected component ofdM. 

Proof. Under changes of the (^-coordinates the last term in Q transforms like a 
Riemannian metric on dM. A change of r-coordinate, i.e. setting r — gi(ip)p + 
g2(<p)p 2 +ff3(</ 3 )p 3 + 0(p 4 ) will give the required (1 + 0(p 2 )) dp 2 behavior iff gi = 1, 
92 = 0, and then r — p + g3(tp)p 3 yields dr = (1 + ip 2 gz) dp + p 3 dg^ and therefore 

(8) (1 - 2r 2 S + 0(r 3 )) dr 2 = (1 - 2p 2 S + 6p 2 g 3 + 0(p 3 )) dp 2 + 2p 3 dpdg 3 + 0(p 6 ) 

so the mixed term remains 0{p i ) (as required in ([7])) iff g 3 is locally constant, and 
the dipidpj term is changed only by 0(p e ), so the metric on dM is well-defined. □ 

The reason for introducing the notion of cusp metric is the following proposition. 

Proposition 2.3. Let X C R" be a space with a cusp singularity of order 2 at 
p. Let g be a smooth Riemannian metric on W 1 and gx its restriction to X \ {p}. 
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Then the pullback g^ of gx to X\ dX , where X is the resolution of X , extends to 
a cusp metric on X . 

Proof. We use the second definition of cusp singularity of order 2, because we can 
use normal coordinates on the first blowing-up space. They do much of the work 
for us. We first work on W 1 and its successive blowings-up and then restrict to X. 
Let p = 0, Z' = [R",0] and Z = [Z',q], with blowing-down maps f3 p : Z' -> E™, 
P q : Z — > Z' , f3 — /3 p o j3 q as in ^ . We need to find a suitable r-coordinate on Z so 
ft* g has the form Q near the interior of the front face. 

Introduce geodesic polar coordinates for g near 0. This means that we choose a 
diffeomorphism of Z' = [M™,0] with R + x S™ -1 so that, with R the coordinate on 
M+, the metric has the form 

P* p g = (dR) 2 + R 2 h(R) 

where R M> h(R) is a smooth family of Riemannian metrics on S™ -1 . In fact, h(Q) 
is the standard metric on S" _1 , but this is inessential for the construction. 

Next let U = (U\, . . . , ?7 n _i) be any coordinates, centered at q, on S n_1 . Then 
h(R) = Ylij h>ij(R,U)dUidUj with smooth. On Z = [Z',q] use projective 
coordinates R,U{ = ^jf-; they are valid in a neighborhood of the interior of the front 
face, so in particular on X near dX. 

Then (i3*h i:j )(R,U) = h V] {R,Ru) = ^(0,0) + O(R) and /3*dUi = d(Rui) = 
Ui dR + R dui and so, with 

(9) S (u) = hjj (0, 0)ujUj , fe 7 = hij(0, 0)dujduj, 

ij ij 

we have 

(10) p*h(R) = S (dR) 2 + RdRdS + R 2 h' + h 

where h = 0(R)(dR) 2 +J2i 0(R 2 ) dRdu,+J2 lJ 0(R 3 ) du,du r Then the main term 

of (3* g = p;p; g is 

(1 + R 2 S) dR 2 + R 3 dRdS + R 4 h! 

and the h part only contributes 0{R 3 )(dR) 2 + 0(R i )dRdu i + Y.ij 0{R 5 )du l du j . 
By the calculation ^ we see that we can make the i? 3 part of the mixed term 
disappear by defining the coordinate r by the relation R = r — \r 3 S, and this 
yields Q, with Cij(r,u) = hij(R, Ru), 

Now consider the submanifold X C Z. Since it is transversal to r = 0, it 
can be locally parametrized as u = u(r, ip) with u v non-singular at r — (here 
if = (ipi, . . . ,ip m -i))- Restricting f3*g to X therefore yields a cusp metric again, 
where S is restricted to X and the dipidipj term is obtained by restricting the 
duiduj-term to X. □ 

Remark 2.4. The proof shows that for a cusp manifold X arising as resolution of 
a space X with cusp singularity p the quantities S and g dx are given as follows. 
First assume that the ambient space is K™ with the standard Euclidean metric. Let 
v be the unit vector at p so that the cusp points in direction v (i.e. the tangent 
cone of X at p is Ryov), and let H be the hyperplane through p perpendicular to 
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v. For small t > let H t = tv + H and X t = -tv |(ln H t ) C H. Then 
dX = lim t _j.o t~ 2 X t C H in the Gromov-Hausdorff sense, and 

(11) S(u) = |u| 2 , = restriction of g eucl to dX 

where g euc i is the Euclidean metric on the subspace H C K™ and |u| is the Euclidean 
norm. For example, if X is given as in |2| then H = R" -1 with the standard metric. 

In general the ambient metric induces the structure of Euclidean vector space 
on the front face of the quadratic (or iterated standard) blowing-up of the ambient 
space, and then S and g a j^ are given in the same way by restriction to dX. 

The dotted lines in Figure [7] indicate the tangent cone of X at p and the line 
u = in its resolution. 

Remark 2.5. The proof shows that the Proposition holds for the more general class 
of metrics obtained from any conic metric by blowing-up of a boundary point. 

3. The rescaled geodesic vector field 

We now investigate the geodesic vector field for a cusp metric Q. We always 
work in coordinates (r, ip) in which ^ holds and write the metric in simplified 
notation as 

g = (1 - 2r 2 S + 0(r 3 )) dr 2 + r 4 2B dr dip + r 4 C dip 2 

where dip = {dipt,. . . ,dip m -i), B = (b%, . . . , 6 m _i), C = (cij)i,j=i t ...,m-i are smooth 
functions of r, <p and S is a smooth function of ip. All vectors will be treated as 
column vectors. By assumption C > at r = 0. Here and in the sequel, 0(r) 
denotes r times a smooth function on DM, hence can be differentiated. 

Let £, r\ denote the cotangent coordinates dual to r, ip, respectively. In order to 
calculate the metric dual to g, we need to invert the coefficient matrix of g. We 
use the general formula (for a scalar a, vector b and symmetric matrix c) 

[b c) = d [-c-H c-H + c-W<ri)> d = a - bc b 

which can be verified by direct calculation, with a = 1 — 2r 2 S + 0(r 3 ), b = r A B, 
c = r 4 C. Then d = 1 - 2r 2 S + 0(r 3 ), hence eT 1 = 1 + 2r 2 S + 0(r 3 ), and we get 

g* = (1 + 2r 2 S + <3(r 3 )) £ 2 + 2 (-B t C^ 1 + 0{r 2 )) £n + r~ 4 (C^ 1 + 0(r 4 )) if 

= e + 2r 2 G(r,<p,Z,*> i ) 

(C^ 1 !] 2 is to be interpreted as rfC -1 ^) where 

G(r, <p, C 9) = G (r, ip, £, 6) + G(r, ip, £, 9), 

U ~' G = S(^e + hC(^- 1 9 2 , G = 0(r)e + 0{r)t9 + 0{r i )9 2 

The energy E(r, ip, rj) is defined as half the dual metric, so if we rewrite it in 
the rescaled coordinates Q we get 

(13) E(r,<p,t,9)= l -e + r 2 G{ Vl r^,9) 

The geodesies are projections of integral curves of the Hamiltonian vector field 
associated with the energy function E and the symplectic form dr A d£ + dtp A dr/, 
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that is they are the r, ip part of the solutions of the following differential equation: 

r' = Ej: £' = -% 
<p' = Er) rf = -E v 



(14) 



In order to rewrite this using the 9 coordinate, we differentiate Q with respect to 
the various variables and obtain 



E 6 

En 



E^ E r — E r + 3v OE^ 
- 3? * -EL 



r^Er, E v 



and 



9' = ^-3^ 



Therefore, in the new coordinates the differential equations (14) are 

a 

r' = E f 



(15) 



r = -E r 



r 



v' ' = \e b e 1 = -Le v -3-e 6 



This new system corresponds to a vector field W that we call the geodesic vector 
field, which is well defined and smooth outside {r — 0}. 

Since E = + r 2 G with G smooth as a function of r, <p, £, 9 by ( 13 1, the rescaled 
geodesic vector field V = rW is smooth up to r = 0, and its integral curves are 
the solutions of the system 



(16) 



r = r[f + r 2 G ? ] £ = r 2 [-2G - rG r + 39G e ] 
ifi = G g 9 = -[Gp+3£9 + 3r 2 6»G 5 ] 



(Throughout the paper the time variables for integral curves of W, V are denoted 
r, t, respectively, and -4- by a prime and 4r by a dot.) Unit speed geodesies lie on 
the energy level {E = |}. This is a (2m— l)-dimensional manifold which intersects 
the boundary {r = 0} transversally and whose intersection with the boundary has 
two connected components corresponding to £ = ±1, by ( |13[ ), In a neighborhood of 
this intersection each component can be parametrized by (r, 99, 0). The f equation 
shows that £ > for geodesies leaving the boundary. We will always consider this 
component. 

It is important to understand the leading terms of the energy and of the rescaled 
geodesic vector field near r = 0. These live on the energy hypersurface {E = 4}, 
which is non-compact since 9 is not restricted to a bounded set (as opposed to rj). 
For the following estimates it is important that they are uniform in all variables in 
a neighborhood of {r = 0}. In particular, 0(r) means r times a smooth bounded 
function on {E = ^}. 

Recall Q. Clearly \S£ 2 + G\ < K{^ 2 + r 2 C- 1 9 2 ), and this easily implies that 

(17) £, 2 + r 2 C~ 1 9 2 = l + 0(r 2 ) on {E = |} 

In particular £ and r9 are bounded on {E — g} and hence also G^, rG r are bounded 
while Gg = 0(r), G v = 0(r) + 0(r9) (the two error terms are not comparable since 
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9 is unbounded) and 39Gg -2G = 2C- 1 6 2 + 0(1) = r~ 2 [2(1 - £ 2 ) + 0(r 2 )] (using 
(17 1 again), so ( fl6| yields, with uniform estimates on {E = |}, 

r = r^ + 0(r 3 ) £ = 2(1 - £ 2 ) + 0(r 2 ) 

P=(G )e + O(r) 9 = -(G ) v - 3£# + 0(r) + O(r0) 

Lemma 3.1. In a neighborhood of {r = 0} i/ie singular locus o/V on i/ie energy 
hypersurface {E =5} is {r = 0, = <S^ = 0}. 




Proof. By ( 18 1 we have £ > 1 — 2£ in a neighborhood of {r = 0}, so £ = implies 
|Cl — ^75 an d hence |f| > -^r + 0(r 3 ), so at a singular point of V we must have 

r = and hence £ = ±1. Then p = C~ x 9 = gives 9 = 0, and then 9 = gives 
5 V = 0. □ 

Therefore singularities of V correspond to critical points of S, that is, to points 
P = P(fo) € dM where S v ((po) = 0. We denote the corresponding singular point 
of V having £ = 1 and 9 = by p e {E = §}. 

Lemma 3.2. T/ie linear part o/V on i/ie energy hypersurface at a singular point 
p = (0,po, 1,0) is 

r = r, ip = 71?' + C _1 0, 6» = 7 2 r - 5^ • (tp - tp ) - 39 

with S vv , C evaluated at p and for suitable 71,72- 

The precise values of 71,72 € M m_1 are inessential for the analysis. 71 = 
~B t C^ 1 and 72 is determined by B,C and the 0(r 3 ) term in the coefficient of 
dr 2 in g. We rewrite the linear part in matrix form: 

/A /i 
(19) Linear part of V: \<fi = 7x 

W \72 S v 

Proof of Theorem |1.2[ The proof is analogous to the proof of the result of 
Melrose & Wunsch in the conic case |MWj . 

If S is constant then every point on C = {(r, <p, £, 9) : r = 0, 9 = 0, £ = 1} is a 
singular point of V. Let p = (0,^0, 1,0) G C. The set {ip = ipo} is transversal to 
C, and projecting the linear part of V at p onto the tangent space of this set we 
obtain 

' 1 
v72 -3 Id, 

(delete the second block of rows and columns in the matrix (19)). This has a 
single eigenvalue 1 and the (m— l)-fold eigenvalue —3, so V is normally hyperbolic 
along C. The unstable eigenvector v (for the eigenvalue 1) has non- vanishing d r - 
componcnt. By the stable manifold theorem [PSj there exists a unique smooth 
unstable submanifold N of dimension m with ON = N n {r = 0} = C, transversal 
to {r = 0}, which is the union of trajectories of V starting at points of C. The 
restriction of V to N is a smooth vector held on N which vanishes at r = 0, hence 
W = 1~V\n extends smoothly from N n {r > 0} to all of TV by Taylor's theorem. 
Also, W has a non-zero 9 r -component at dN by the f equation in (15). Near 
r = 0, N is a graph over (r, ip) since its tangent space is spanned by TC and v. 
Therefore, the projection (r, ip, 9) 1— > (r, ip) restricts, near dN, to a diffeomorphism 
ip : N -> M. Therefore, the vector held W = ip*(W\ N ) is well-defined near dM, 
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and its integral curves are geodesies in M. Since W is transversal to dM, the 
result follows from standard facts about flows. 

Remark 3.3. l)This can be generalized to the case when S is constant only on an 
open subset of dM , with the same proof. 

2) By the standard Gauss lemma (which is also applicable here, cf. the footnote in 
the proof of the conic theorem in |MW| ), the theorem implies that for constant S 
the normal form Q may be improved: There is a trivialization of a neighborhood 
U of the boundary, dM x [0,To) = U (given by the flow of the vector field W in 
the proof), so that 

g {u = dr 2 + r A h{r) 
where h is a smooth family of metrics on dM . 



4. Local behavior of the rescaled geodesic vector field near 

singularities 

Let p be a singular point of the rescaled geodesic vector field V corresponding 
to a critical point p of S. We will determine the local phase portrait of V near p 
in the case where p is a non-degenerate critical point of S. 



We first determine the eigenvalues and eigenvectors of the linear part ( 19 1 of V 
at p. Recall that the Hessian of S at a critical point p is a well-defined symmetric 
bilinear form on T p dM. By way of the metric gdM, this form corresponds to a 
linear endomorphism of T p dM , which we call the Hessian of S relative to gaM at 
p. For example, if one chooses coordinates in which d Vl , . . . , d Vm _ 1 are orthonormal 
with respect to gdM at p then this is simply the matrix of partial derivatives S (fii , Pj , 



and C = Id at p. A simple calculation using ( 19 ) yields 



Lemma 4.1. The linear part o/V at a singular point p has an eigenvalue A^p) = 1 
with eigenvector transverse to r — 0, and for each eigenvalue a of the Hessian of 
S relative to gdM a t P an d each corresponding eigenvector u-J^ — Ui T^pl 

the 

two eigenvalues 

(20) A 2 (p,a) = and A 3 (p,a) = 

with eigenvectors 
(21) 

d d d d 

i>2(p,a,u) = u— + \ 2 (p,a)(Cu)— and v 3 (p,a,u) = u— + X 3 (p,a)(Cu) — 

Note that the eigenvalues A2 ; 3(<^, a) are non-real if a > 9/4. 

This allows us to analyze the geodesies near non-degenerate critical points. We 
only note the cases of maxima and minima of S, other critical points can be treated 
similarly. 

Proposition 4.2. Let p be a non-degenerate critical point of S on a cusp manifold 
M. 

(1) If S has a local maximum at p then there is a neighborhood U of p in M 
so that U H M is foliated by geodesies. 

(2) // S has a local minimum at p then there is a unique geodesic which starts 
at p. 
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Proof. At a local maximum of S, all eigenvalues a of its Hessian are negative, 
so all eigenvalues of the linear part of V are real and A3(p,a) < < Aa(p, a). 
By the unstable manifold theorem, the unstable manifold M£ is a smooth Tri- 
dimensional submanifold of 3 T*M whose tangent space at p is spanned by the 
eigenspaces corresponding to Ai and the various A2(p, a). Eigenvectors for Ai have 
non-zero <9 r -component, and the (^-components of the eigenvectors 1*2 (p, a) span 
T p dM . Therefore the projection of M™ onto M is a diffeomorphism in a neighbor- 
hood of p. Since Mp \ p is locally foliated by trajectories of V, the result for a 
maximum follows. 

At a minimum of S, all values a are positive, so the linear part of V has a unique 
eigenvalue, Ai = 1, with positive real part. By the unstable manifold theorem, 
there is a unique trajectory leaving p. The result follows. □ 

From now on we focus on cusp surfaces. In this case there is only one value of 
a, which we denote by 

o(p) - ^(P)C(P)- 1 
In arc length coordinates this is simply S vv (p). We write A2,3(p) = A2,3(p, a(p)) 
and 1/2,3 (p) = 1/2,3 (p,a(p),l). 

Remark 4.3. In the case where M = X is the resolution of a surface with cusp 
singularity, the geometric meaning of a(p) is as follows: Recall that dM C R™ -1 , 
where K ra_1 is equipped with the standard euclidean scalar product, S(u) = |u| 2 and 
g9M * s the induced metric. Then 

§a(p) = 1 + p ■ tf(p) 

where K(p) is the curvature vector of dM at p. This is easily seen by taking arc 
length parametrization ip H> u((p), then S vv — 2(\u ip \ 2 + u • U w ) and \u v \ = 1, 
u w = K . Now evaluate at u = p. 

Linearization near a singular point. A triple of pairwise distinct complex num- 
bers Ai,A2,A3 are resonant if there exist i £ {1,2,3} and 01,02,03 £ No with 
a,i + 02 + 03 > 2 such that 

(22) aiAi + a 2 A 2 + 03A 3 = A^ 

This relation, or the triple (ai, 02, 03), is called a resonance relation of type i among 
Ai, A 2 and A3. 

Given a non degenerate local extremum p of the function S, we observe that we 
always have the relation 3Ai(p) + A 2 (p) + As(p) = 0. Multiplying by any natural 
number and adding Aj(p) we obtain three families of resonance relations. 

We recall the following very useful linearization criterion (Samovol's criterion). 

Theorem 4.4 f |Sa|). Let (/Xi . . . (i s , o~\, . . . , a u ) be the set of eigenvalues of the 
linear part of a smooth vector field V at 0, a singular point of V . Assume that 

Re(ps) < < Re{m) < < Re(ai) < ...< Re{a u ) 

Assume there exists a positive integer I such that for each resonant relation, 

X = E^/^ + E^jj 
where A is an eigen-value, there exists an integer i < s or j < u such that 
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—lRe{ni) < —Re(r^ii\ + . . . + r^Mi) or lRe(cj) < Re{r^a\ + . . . + rfo-j). 

Then the vector fields V is C -conjugate to its linear part in a neighborhood ofO. 

Remark 4.5. The condition A 2 (p) 7^ —1 is equivalent to a(p) =/= 2. 

As a consequence of this result and of the remark, we get the following proposi- 
tion. 

Proposition 4.6. Assume that M is a cusp surface and that p £ dM is a non- 
degenerate local minimum for the function S for which a(p) < | and a(p) 7^ 2. // 
•^(p) € Q then V is C 1 linearizable near p, otherwise it is C 2 linearizable. 

Proof. Write Ai — Xi(p). The assumptions imply Ai = 1, A3 < — | < A2 < and 
A 2 ^-l. 

1) Assume A2 is irrational. The resonance relations between the eigenvalues are 
of the form 6(3, 1, 1) + R, where R £ {(1, 0, 0), (0, 1, 0), (0, 0, 1)} and 5 £ N. For 
each such relation we always find 

2Ai < [36 + w]Ai for ui = or 1. 



We conclude using Theorem 4.4 that V is C linearizable at p 



2) Assume that A2 = — | 7^ — 1 is a negative rational number written in its 
irreducible form. 

Let (ai, a 2 , 0,3) be a resonance relation of the form aiAi + 02X2 + CJ3A3 = \%. 
Then a\ > 1 since otherwise > CI2A2 + (I3A3 = (1 — ai)Ai > 0. Hence aiAi > Ai. 

So for all resonance relations of type 1 we always find aiAi > Ai, so that we can 
take j = 1, 1 = 1 in regards of Samovol's criterion. 

Let (ai, 0,2, 03) be a resonance relation of the form aiX\ + a 2 A2 + a 3 A3 = A 2 . 
If ci2 > 2 then a2(— A 2 ) > — A2. 

If ci2 = 1 then ai + a3A3 = 0, thus a\ > 2 and so aiAi > Ai. 

If 02 = we check that (01,03) 7^ (1,1) since A 2 7^ —1. Thus we either have 
aiAi > Ai or a 2 (-A 2 ) + a 3 (— A3) > -A3. 

So for all resonance relations of type 2 we can take I — 1 and, respectively, i = 1, 
j = 1, j = 1 or i = 2, in view of Samovol's criterion. 

Let (a±, 02, a 3 ) be a resonance relation of the form a^Ai + a 2 A2 + asA3 = A 3 . 
If aa, 03 > 1 then 02(— A 2 ) + 03(— A3) > —A3. Note also that the cases a 2 = 03 = 
and et2 = 0, 03 = 1 are impossible. 

If 02 = 1, 03 = then A3 = aiAi + A2 > A2, which is impossible by hypothesis. 

So for all resonance relations of type 3 we can take i = 2, 1 = 1, in view of 
Samovol's criterion. 



We conclude using Theorem 4.4 that V is C 1 linearizable at p. □ 

5. Global boundary dynamics 

The rescaled geodesic vector field V is tangent to the boundary r = of {E = 
5} C 3 T*M, which is the union of two components corresponding to £ = ±1. We 
will study the dynamics of the flow of V on the part £ = 1. As explained in the 
introduction this can be identified with the space T*dM with variables ip, 9 in local 
coordinates. 
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Let 

E d &, 9) = S(<p) + \C-H 2 = (G )| C=1 . 

By (18) the trajectories of the vector field V on T*dM are the solutions of the 
system 

(23) <p = El, B = -E°- W 

Explicitly, 

y = CT% = -S v (<p)-2>6-\{C- l ) v 9 2 , 



The system (23) is a Hamiltonian system modified by the damping term —39. It 
describes motion of a particle on DM, with Riemannian metric given by ggjw, m the 
potential S and in the presence of damping (for example, friction) . The following 
properties are standard for such systems. We provide a proof for completeness. In 
this discussion the dimension of M is arbitrary. 

Proposition 5.1. Let M be a cusp manifold. The singular points o/V on T*dM 
are equal to the critical points of the boundary energy E d and given by the critical 
points of S: 

{(<p,9): 9 = 0, S v ^)=0}. 

Every maximal trajectory 7 o/V on T*dM is defined for all times. If S is constant 
or has only isolated critical points then we have in addition: 

a) As t — > 00, the trajectory "/(t) converges to a singular point o/V. 

b) As t — > —00, either E 9 is bounded along 7 and j(t) converges to a singular 
point o/V, or else E 9 (^/(t)) — > 00. 

Proof. Let t j(t) = ((p(t),9(t)), t e (T_,T + ) be a maximal trajectory of V. 
From ( 23 ) we get 

(24) E 9 = -Z\9\ 2 

where \9\ 2 := C~ x 9 2 . In particular, E d is decreasing along 7. Since E d is bounded 
below, it must approach a limit as t — > T + . In particular, 7 is contained in a 
compact subset of T*dM, which implies T + = 00. Fix T £ (T_,oo). Equation 
(24) implies \9{t)\ 2 dt < 00, and then the boundedness of 9 easily implies that 
9(t) — >• as t —> 00. Let a(t) = S v ((p(t)). The equation for 9 and 9 — > imply that 
a(s) ds — > as t — > 00, and boundedness of a again implies that a(t) 0, so 
7 converges to the set of critical points of E d . If this set is discrete then a) follows 
immediately. If S is constant then (C -1 )' = (^(C -1 )^ = 0(9), hence 

{C- 1 9y + '6C- 1 9 = 0{\9\ 2 ) + G- x 9 + 3C- X 9 = 0{\9\ 2 ) - iCT^Cr 1 )^ 2 = 0{\9\ 2 ) 

which is integrable over [T, 00). Now 9 — > implies that limT'-^oo fx {C~ 1 9)' dt 

exists, hence so does limr'^oo 3C~ 1 9 dt and hence lim^oo tp(t) exists. Choosing 
T sufficiently big one may assume that all these curves lie in a single coordinate 
patch. This proves a). 

Next, §24§ implies E 9 = -6E 9 + 0(1) which implies \E 9 (-f(t))\ < ife 6|t| for 
some constant K. Then the same estimate must hold for |#(i)| 2 , and this implies 
T_ = —00. Now by monotonicity either E 9 (^(t)) — > 00 as t — > —00, or else it is 
bounded. In the latter case (24) implies that f_ \9(t)\ 2 dt < 00. Then the same 
arguments as in the proof of a) yield b). □ 
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Figure 3. Boundary dynamics between a maximum of S (at </? = 
0) and a minimum in case T = oo: The dotted line is the graph of 
the barrier function /, the solid line is the trajectory T p 



Remark 5.2. The proof carries over almost literally to the following perturbed 
version of b): Let 7 be a solution of the system 

tp = E$ + f{t), 9 = -E°-39-g(t) 

where \ f(t)\ + \g(t)\ < Le^ /K for all t < 0, for some constants K, L. Then, if E 9 
is bounded along 7, then j(t) converges to a singular point o/V as t — » —00. 

We now analyze the case of surfaces M in greater detail. Then we may and will, 
for simplicity, assume: 

(25) tp is an arc length parameter for dM. 

This means C = 1 on dM, so V = 9d ip — (S^ + W)8q when restricted to dM. 
For the following discussion we assume: 

S has a maximum at ip — with S vv (0) < 
^ ^ and a minimum at p m i n > 0, and S v < on (0, <^ m i n ). 

Denote by p-,p m i n the points in dM given by p — 0, tp — tp m i n , respectively. 
Recall that these correspond to singular points p_,p min of V in T*dM, given by 
these values of tp, and 6 = 0. Recall that at a local maximum of S, the linear part of 
V restricted to T*dM has eigenvalues A3 < < A 2 , with eigenspace for A2 spanned 
by vi = d v + A2<9e- Thus, by the unstable manifold theorem there are (up to time 
shift) precisely two trajectories starting at p_ at t = —00, and their directions at 
p_ are ±1^. We discuss only the trajectory starting in direction v-i and denote it 
by T p , there is an analogous discussion for the other trajectory. 

The trajectory 1 1-> T p (t) = (tp(t), 9(t)) starts out at p_ into the first quadrant 



tp > 0, 6 > 0. It is defined for all times by Proposition 5.1 a). As long as tp < <p m i n , 
that is, for t < T := sup{£ : tp(s) < tp m i n for all s < t}, it will stay in the first 
quadrant since V points inside this quadrant along this part of its boundary, more 
precisely: 

V(tp, 0) = -S v (tp)d 6 and - S 9 (<p) > for < ip < tp min 
V(0, 6) = 6d v - Wd e and 6 > for the first quadrant. 



See Figure [3] Since tp = 9 > in the first quadrant, tp is strictly increasing on 
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(— oo,T). Now there are two possibilities: 

either T < oo, then tp(T) = tp m[n and 8{T) > 

or T = oo, then ip(t) < p m [ n Vi and ip(t) — > (f m i n , 0(t) — > as t — > oo 



(and hence T p (t) — > p m ; n ), by Proposition 5.1 a) since there are no critical points 
with ip G (O^nun). 

We discuss the first case below in the proof of Theorem |1.5[ In the second case 



in (27) the trajectory T p is the graph of a smooth function h : (0,p m i n ) — ► (0, oo). 
We now give a sufficient condition for this to occur. It is obtained by a simple 
barrier argument. 



Lemma 5.3. Suppose we are in the setup described above, starting from (25 1. 
Assume that there is a C 1 function f : [0, ip m i n ] —> M+ satisfying the conditions 

(28) f(v>min) = 0, and f((f) > for tp € (0, cp min ) 

(29) /(0) > 0, or /(0) = and /'(0) > A 2 (p_). 

(30) //' + 3/ + S , „>0on(0,¥w) 
T/ien: 

a) TTie closure of the image ofT p _ is £/ie graph of a C 1 function h : [0, <p m in\ ~> 
E+, and h'(0) = A 2 (p-). 

b) If in addition (//' + 3/ + S 1 ^)' < ai <p m m i/ien T p _ approaches p m in 
tangent to the eigenvector ^(Pmm), that is, h'(ip m i n ) — A2(p m m)- 

Concerning the condition in b), observe that f(<f m in) = implies differentiability 
of //' there, with derivative (/'(<£ m in)) 2 - Also, the function g = ff + 3/ + S v 
satisfies g(ip m i n ) — and g(ip) > for tp < tp m i u , which implies that <?'(<p m in) < 0. 
The given condition strengthens this to g'(<p m i n ) < 0. 

Proof. Let 

D(f) = {(<p,0):Q<<p< cp min , 0<e< f(cp)} 
We claim that the vector field V = 6d v — (S v + 39)dg points inside D(f), or is 
tangential to the boundary, at each boundary point of D(f): This was checked 
before the proposition for the parts of the boundary where 9 = or ip = 0. The 
remaining part of the boundary is the graph of /, which is the zero set of f(ip, 0) = 
f(<p) - 6, with / > inside D(f). Now V/ = Of + S + 30, which at = f(tp) 



equals //' + 3f + S, so (30) gives V/ > which was to be shown. 

By (29) the trajectory T p _ starts out into the interior of D{f), and by what we 
just proved it can never leave it. Since ip = is positive in the interior of D(f), 
the trajectory traces out the graph of a smooth function ip i-> h(tp) on (0, </? m in)- h 
extends to be smooth at by the unstable manifold theorem, and it is C 1 at ip m [ n 
since the trajectory is tangent to an eigenspace at p min . This proves a). 

Finally, (//' + 3f + S v )' = (/') 2 + 3/' + S lfnf> at ip m i n , and if this is negative then 
/' (everything evaluated at >p m i n ) must lie strictly between the zeroes A2 and A3 
of the characteristic polynomial A 2 + 3 A + S vv> of the linear part of V at (p m i n , so 
h! > f > A3. Since h! must be either A2 or A3 by the local behavior of trajectories, 
it follows that h' = X 2 and hence b). □ 

Proposition 5.4. Let M be a cusp surface and assume S is a Morse function on 
dM satisfying the bound S vip < 9/4 everywhere, where p is an arc length parameter 
on dM. Let T — (J T p where p ranges over the critical points of S and T p is the 
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unstable manifold of the point p for the vector field V on T*dM. Then T is the 
graph of a C 1 section of T*dM . That is, 

T = {(ip,9) : (p € dM, 6 = h((p)}, h e C\dM). 

Furthermore, h'(ip) — \2(p(<p)) at each critical point p = p(ip) of S. 

In fact, the proof shows that h is C°° except at minima of S, and at the minima 
the regularity can, in general, be slightly improved, see the remark after Proposition 



Proof. Since S is a Morse function and dim dM = 1, its critical points are either 
maxima or minima, and these alternate along dM. We first consider an interval 
between a maximum p and the next minimum p m i n . We may assume that the 
maximum is at ip = and the minimum is at ip m - ln > 0. Then we are in the 
situation of Lemma p. 3| As barrier function we take / = — This clearly 



satisfies condition i28f. Also, condition (301 holds since 



Q '• — f f 3f -\- Sip — ^S^S^p^p -\- — Stp^gStptp 1) > 



on (0, (p m in), and .g'l^min) < since S v = 0, - 1 < 0, S vv > at <p n 

Finally, since /(0) — condition (29 1 is equivalent to — |a > ~ 3 +^ 9 ~ 4a . where 



a = S vlp (0), which is easily checked to be true for any a < 0. So Lemma 5.3 a) and 
b) are applicable and yield a function h p _ Pmiii on the interval [0, ^ m i n ] satisfying the 
claim of the Proposition there. We construct corresponding functions for each pair 
(p_, p m in) of consecutive maxima and minima of S, and by changing the signs also 
for pairs (p m i n ,p+) of consecutive minima and maxima; here ft. PminP+ is negative 
in the interior of this interval. All these functions fit together to form a continuous 
function h on all of dM since they vanish at the endpoints of these intervals, h 
is smooth at maxima p of S since its graph is the unstable manifold of p near p, 
which is smooth, h is C 1 at minima p m i n of S since its derivatives from both sides 



agree, by Lemma |5.3p ). □ 

The value | is optimal for this type of barrier function. Note that the condition 
< S vl p < | at a minimum p m i n is equivalent to the linear part of V at p m i n 
having two distinct negative real eigenvalues. 

We now consider the case where M = dX is the resolution of a cusp singularity. 
Recall from Remark 2.4 that dX is naturally the subset of a Euclidean vector space. 



If X is given as in (|2j) then this is simply K" 1 x {0} with the standard metric. In 
general we still may identify it with K n_1 with the standard metric. 

Proposition 5.5. Let X be a surface with cusp singularity. Assume that dX C 
M n_1 is contained in the boundary of a strictly convex set which contains the origin. 
Also, assume that, for any point p € dX where dX is tangent to the sphere through 
p centered at the origin, it is only simply tangent to that sphere. Then we have 
Stp V < 2 everywhere on dX , where tp is an arc length parameter. In particular, 



Proposition 5.4 and Theorem 1.4 apply to M = dX 



In the case of a surface in M 3 the simple tangency condition is equivalent to the 
condition that the osculating circle of dX, wherever it is defined, be never centered 
at the origin. Strict convexity is meant in the sense of nonzero curvature. 
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Proof. Choose an arc length parametrization p H> u(<p) of dX. Then S(<p) — 
\u(ip)\ 2 , so the simple tangency condition is equivalent to S being a Morse function. 



By Remark 4.3 we have \S vv (p>) = 1 + u • K(u) for all u = u(tp) where K(u) is 
the curvature vector of dX at u. Suppose dX is contained in the boundary of the 
strictly convex set ]C, and let P u be a supporting hyperplane for K. through u. Then 
K(u) points into that closed half-space determined by P u which contains JC. Since 
JC contains the origin, it follows that u- K (u) < 0, so S vip (p) < 2 < | for all ip. □ 

The proof shows that the conclusion also holds with a certain amount of non- 
convexity. 



6. The exponential map 

In this section we prove Theorem define the exponential map and then prove 
Theorems Ol Ol and [L5l 



Proof of Theorem The geodesies are, up to reparametrization, the projec- 
tions to M of the integral curves of the vector field V on the energy hypersurface 
{E = |} C 3 T*M. Thus, let 7 : (T'_,T' + ) -> 3 T* M be a maximal integral curve 
of V, and assume that 71/7^ t') ^ s contained in {r < ro} for some T' € (T!_,TV). 
The number rg will be chosen sufficiently small below. Then necessarily T'_ = —00 
since V is tangent to the boundary. We will show that lim t _ y _ 00 7(f) exists and is 
a singular point of V. 

Write j(t) = (r(t),(p(t),£(t),0(i)). Here we have chosen a trivialization of M 
near the boundary, which allows us to identify 3 T*M with T*R + x T*dM there, and 
then (r, £) denotes a point in the first factor and (<p, 9) a point in the second factor. 



The proof falls into two parts: First, we use the f and £ equations in (18 1 to show 
that £(t) must be close to one on (— 00, T"), hence r(t) exponentially decreasing 
as t decreases, and that 6(t) is bounded. Then the decay of r together with the 
boundedness of 9 shows that the ip, 9 dynamics is close to the boundary dynamics, 
which allows us to use Proposition |5. lb ) 



By ( 18 1 there is K > such that £ > 2(1 — £ ) — Kr on the energy hypersurface 



Claim 1: £ > \ for t < T . 

Proof. For ro sufficiently small we have £ > 1 for |£| < |, so if £ (to) < \ for some 
t < T' then £(t) < -\ for all t < to — 1, and then f = £r + 0(r 3 ) implies that 
r is unbounded as t — > —00, contradicting the assumption that r(t) < ro for all 
t < T'. □ 



Claim 2: £ 2 > 1 - Kr 2 for t < T'. 

Proof. Let f(r,£) = £ 2 + Kr 2 . The part of 7|(_oo,T') where / > 1 is forward 
invariant under the flow of V since at any point with f(r, £) = 1 we have 2(1 — £ 2 ) — 
Kr 2 = Kr 2 and hence |V/ = &+Krr > ZKr 2 + Kr(r£+0(r 3 )) > Kr 2 + 0(r 4 ) > 
0, by Claim 1, for ro sufficiently small. Hence if / < 1 at some point 7(^0) then this 
remains true for all t < t . Now /(r,£) < 1 implies 2(1 - £ 2 ) - Kr 2 > 1 - £ 2 > 0, 
so we would get £ > 1 — £ 2 > for all t < t . Then £ would be monotone 
increasing in t < to, so \im t ^_ OD ^(t) would exist and lie in [|, 1). However, this is 
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impossible since it would imply £(io) — 00 ) > /_° (1 — £,(t) 2 )dt = 00 and hence 
unboundedness of £. □ 

From Claim 1 and r = r£ + 0(r 3 ) we obtain r(t) < ife~'*'/ 3 , and Claim 2 gives 



boundedness of 6{t), using ((TtJ) . Also, ((T7| implies that - 1| < K' e" 2 ^^ . 
Therefore Remark 



5.2 



applied to the p and 9 equations in (18), implies that j(t) 



approaches a singular point of V as t — > — 00. 

The claim that each singular point p of V is the starting point of a trajectory of 
V follows from the invariant manifold theorem and Lemma 14 . 1 1 since the unstable 
tangent space at p contains the eigenvector for the eigenvalue A = 1, which points 
inside M. 

Finally, returning to the unrescaled geodesic vector field we have r' = £ + 0(r 2 ) 
along a geodesic, and since £ is near one this implies that the geodesic reaches r = 
in finite backward time, so T_ is finite. This completes the proof of Theorem |1.1| 

Definition of the exponential map. We now define the exponential map in the 
case where M is a surface with connected boundary and S is a Morse function on 
dM. Then dM is diffcomorphic to a circle and the critical points can be labeled, in 
cyclic order, Pi, pi, P2, P21 ■ • • j Pfc, Pk f° r some k € N where S has maxima at the p.; 
and minima at the p^. This fixes an orientation of dM. Recall the considerations 



in the proof of Proposition |4. 2| in the case of surfaces: The unstable manifold M" is 



two-dimensional for each maximum p = Pi, and it projects diffeomorphically to M 
in a neighborhood of p. It is one-dimensional for each minimum p = p^. Then for 
each i the integral curves of V starting at p i: considered as sets, foliate M p . \ 
They may be parametrized by the points of a small closed half circle Ti around p, ; 
in Mp. transversal to V. Identify Tj with the closed interval [i — 1, i], so that the 
parametrization preserves orientation and the boundary points correspond 
to curves in the boundary. In addition M™, \ p^ consists of a single interior integral 
curve, which we label by i. In this way every interior integral curve starting at a 
boundary point is labeled by a unique number in (0, k]. Identifying S 1 with K/(fcZ) 
we obtain a parametrization of the set of geodesies starting at the boundary by S . 

Definition 6.1. The exponential map is defined by the parametrization introduced 
above, together with arc length parametrization along each geodesic starting at t — 
0. 

Note that the exponential map preserves cyclic ordering, with respect to the 
chosen orientation of dM and the standard orientation of S 1 , in the following sense. 
If a, b, c G S 1 are pairwise different and lie in this order on S 1 then there is tq > 
so that the geodesies 7 Q ,7b,7 c do not intersect and lie in this order. Here, To may 
depend on a, 6, c, and for certain cusp metrics will in fact not be bounded away 
from zero, which leads to discontinuity of the exponential map, see Remark 1 6. 5 1 



Remark 6.2. The domain of the parametrization q > j q may be described some- 
what more naturally with the help of the following construction: For a vector field 
V on a smooth manifold S with an unstable critical point p there is a natural no- 
tion of blowing-up of p in S with respect to V. This is a smooth manifold with 
boundary, denoted [S,p]v, together with a blowing-down map /3 P : [S,p]v — > S. 
The boundary (front face) parametrizes the integral curves o/V starting at p. It is 
diffeomorphic to a sphere and may be thought of as small sphere around p transver- 
sal to V. The map f3 p may not be smooth, but V lifts to a smooth vector field on 
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[S, p]v> given by sd s near the boundary for a suitable boundary defining function 
s. The construction generalizes to the case where £ is a manifold with boundary, 
p G <9E, and V is tangent to <9£ ; and then yields a manifold with corners. See 
[HMV1 Section 2] for details. 

We apply this to the surface M™ ( for each i and obtain [Mp. , pjv- Then the 
transversal T, is naturally the front face /3p (p,j) of the blowing-up. The various Ti 
are now glued at their endpoints, which are also identified with the points pj, and 
this is the first factor S 1 of the domain of exp dM . 



Proof of Theorem |1.3[ We constructed the exponential map above. We need to 
show that it is surjective. For this we prove that the union of the unstable manifolds 
Mp over all critical points projects onto a neighborhood of dM under the canonical 
projection tt : 3 T*M -> M. 

First, observe that each unstable manifold M" intersects the boundary {r = 0} 
transversally near p since this is true for the linear part of V, and since both Mp 
and {r = 0} are invariant under the flow, the intersection is transversal everywhere. 
The intersection Mp n {r = 0} is the unstable manifold Tp of p of V restricted to 
the boundary. 

The projection of the unstable manifolds M", are (closures of) single geodesies 
starting at p£, which divide a neighborhood of dM into connected components Ui, 
each containing one maximum p^. It suffices to show that 7r(Mp ( ) contains Ui for 
each i. 

For this, first consider the intersections with the boundary. Fix i and write 



p = Pi and p min = p-. By the discussion before Lemma 5.3 the part T p _ of T p 
leaving p in the positive ip direction projects onto the (open) interval from p 
to Pmin. A similar statement holds for the interval from p' i _ 1 to p_. Therefore, 
H {r = 0} projects onto (7, n {r = 0}. Since Mp. intersects the boundary 
transversally, it projects onto some neighborhood of the open interval from p[_i to 
p'i- We need to show that this neighborhood cannot shrink to zero width when one 
approaches p' i _ 1 or p^. We will prove that Ui intersected with a neighborhood of 
p- is contained in 7r(Mp.); the argument at p' i _ 1 is then analogous. 



Recall the dichotomy (27). If T < oo here then 7r(r p _) contains p[, so we are 
done. On the other hand, if T — oo then r p _(f) — > p min as t — > oo, and the 
eigenvalues A2(p m in), A3(p m j n ) are real. In this case, the behavior of M^_ near 



p min may be understood by linearizing V near p min . We use Lemma 6.3 below, 
applied as explained after its statement. For later purposes the lemma is stated 
more strongly than needed here. Here we only need the consequence that, in a 
neighborhood of p min , the curve Mp m . n is contained in the closure of Mp_, so the 
projection of Af" n is contained in the closure of Ui. This completes the proof of 
Theorem 11.31 



Lemma 6.3. Let A, /i < and let V be the linear vector field x\d Xl + Xx2d X2 + 
HX3d X3 on ]R 3 . Let m > 1. Let S be an invariant C m surface with boundary 
r = S n {x\ = 0}, and assume this intersection is transversal and T is a single 
trajectory of V . Denote by X% the non-negative x\-axis. 

Then a neighborhood £' of zero in SUXi is a surface with corner. More precisely, 
let t = — y and p = ^, so r, p > 0. Assume T C {x2 > 0}. Then £' is a graph 

(31) S' = {(xi,x 2 ,x 3 ) : x 3 = A(xixl)x%, (xi,x 2 )eU'} 
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where A is a C m function on [0, e) for some e > 0, and U' is a neighborhood of the 
origin in the quarter plane M + x K + . 

The lemma is applied as follows, in the context of the discussion before Lemma 



5.3 Let p m in be a minimum of S, and suppose the second alternative in (27) holds. 
Then T p (t) approaches p min as t — » oo. By Proposition |4.6| we may linearize 
V near p min by a C m -diffeomorphism where m € {1,2}. Thus, we introduce 
local coordinates Xi,X2,Xs so that p min is the origin and V is given by x\d Xl + 
Xx 2 d X2 + fixsd x3 on a neighborhood U where {A, ^} = {A 2 (p m m), A 3 (p m i„)} and 
the boundary is given by {xi = 0}. We take T = T p and S = M p _. Clearly 
Afp m . n = X\. Since T is a trajectory of V approaching the origin, it must have a 
tangent vector there, and we may assume that the coordinates are chosen so that 
this vector is d X2 . Then L C {x 2 > 0} near the origin, so the assumptions of the 
lemma are satisfied. 

Proof of Lemma [Q| Choose a point p of T. Since the intersection of E with {xi = 
0} is transversal we may choose a C m -curve u(s) = (s,u)2(s),u>3(s)), s6l = [0, e) 
contained in E, with w(0) = p, so w 2 (0) > 0. Also, since x\x\~ is constant along any 
integral curve of V, only trajectories passing through this curve will contribute to 
E' if this is chosen sufficiently small. 

For s € I let M + —> M 3 , M T s (t) be the forward integral curve of V starting at 
w(s), so r s (0) = u)(s) and F = F, up to time shift. 

For each s, the quantities x\x\ and x-iX~^ p are constant along r s , that is, for 
points (xi,X2,xa) in the image of T s . Evaluating at t — shows that, along T s , 

(32) xixl = soj 2 (s) t , X 3 X2 P — uj 3 (s)uj2{s)~ p ■ 

Now w 2 (0) > 0, hence the equation a — sw 2 (s)' r can be solved C m -smoothly for 
s = s(cr), for a in a half neighborhood J of zero, and the function A : J — > M, 
a i — ^ W3(s(cr))a; 2 (s(cr)) _p is C m . Then we have along r s , with this cr, 

£3 = A(a)x p 2 = A(x\x T 2 )x p 2 . 

The claim follows. □ 



Proof of Theorem |1.4[ The heart of the proof is that the various unstable man- 
ifolds fit together nicely. So we first prove the following proposition. 



Proposition 6.4. Assume the setting and the conditions of Theorem \l.J\ Let M u 

be the union of the unstable manifolds M™, where p ranges over the critical points 
of S. Then M u is a C 1 manifold, and there is ro > so that the intersection of 
M u with {r < ro} projects diffeomorphically to M (1 {r < ro} under the projection 
7T : 3 T*M -> M. 

See Figure |4j 

The proof actually gives more regularity than C 1 : The manifold M u is C°° 
except at points on the curves M™ where p is a minimum of S, and here the 
regularity is C°^ p ' where 

(33) 

with A 2 = A 2 (p), A 3 = 
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Figure 4. The dynamics near the boundary. 



Proof of Proposition \6^4\ The main task is to prove that the different fit to- 
gether to form a C 1 manifold, the projection statement will then be seen to be an 
easy consequence. 

Since integral curves do not intersect, the manifolds M£ do not intersect for 
different p. Since M u is invariant under the flow, it is enough to prove the regularity 
near r = 0. 

First consider the boundaries = dM£ — M£ H {r = 0} and their union 
r n = UpTp- These are the unstable manifolds of V restricted to the boundary 
T*dM . By assumption S w < | everywhere, so we may apply Proposition 



5.4 



and 

conclude that T is the graph of a C 1 function T = {{ip, h(ip)) : cp g dM}. Thus, 
both the regularity and the projection statement of Proposition |6.4| are true for the 
boundary of M u . 

We now prove that M u is a C 1 manifold in a neighborhood of r = 0, with 
boundary T u . Since is a smooth surface for each maximum p of S, we only 
need to prove the regularity in a small neighborhood U of p min where p m i n is 
a minimum of S. Thus, let p m i n be a minimum of S and let p , p + be the 
maxima of S closest to p m i n , so that the corresponding parameter values satisfy 
<p- < Vmin < <P+- Then the set M u H U is the union of M£_ n U, M£ mln n U 



and Afp + fl U. We now apply Lemma |6.3|as explained after its statement, first to 



S_ = Mp_ n U and then to S + = M£ n U. By the last statement of Proposition 
BTil T p± have tangents ±^ 2 (p m in) at p min , so we take A = A 2 (p min ), (i = A 3 (p min ) 



in Lemma 6.3 and this implies p > f there. The lemma shows that H± arc graphs 
as in (31), with possibly different functions A±. 



The more precise regularity statement ( 33 ) follows from the regularity in Lemma 



6.3 together with the regularity of the linearization, Proposition 4.6 

Finally we prove that M u projects diffeomorphically to M near r = 0. For any 
boundary point q g T u = dM u , the tangent space T^M U is transversal to {r = 0}, 
since this is already true for each M£ separately. In addition, TqM u D TqT", 
so TqM" fl Tq{r = 0} contains a vector with nonzero <9 v -component. These two 
facts combine to show that, for each q 6 dM u , the tangent space TqM u projects 
isomorphically to T q M under dn, where ir : 3 T*M — ¥ M is the projection, and this 
implies the claim. □ 

We now prove Theorem [Oj By Theorem [Tl] every geodesic starting at dM 
must do so at a critical point p of S, so the corresponding integral curve of V 
starts at a singular point p of V. This integral curve is then contained in the 
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unstable manifold Mp of p. Conversely, n {r > 0} is the union of such integral 
curves. Therefore, Proposition 6.4 implies that exp dM is bijective. 

It remains to prove the continuity of exp aM and of its inverse. We continue to 
work on M u . Some care needs to be taken since the vector field W blows- up at 
the boundary. But the main point is that only its d v component blows-up, while 
its d r component is £ + 0(r 2 ) = 1 + 0(r 2 ). Recall that the domain of exp aM is the 
set of (q, t) where r G (0, To) and q is a point on a transversal near p, for some 
i, where the transversals are glued at their endpoints. The continuity of exp 9M at 
points (qa,r) where go is not a boundary point of a transversal 7* is clear, so we 
assume that qo is a boundary point of Tj, say the 'right' boundary point, which 
labels the geodesic starting at the next minimum p^. The idea is this: As sets, the 
geodesies j q , with q in the interior of 2j approaching qo, converge to the union of 
the boundary trajectory r p; and of M™, . Since r' = ^ is approximately one, the 
part of 7 9 near T Pi is traversed in a very short time. Then the part near M™, must, 
including its time parametrization, be close to 7 9o . 

More precisely, fix (qo, t) and a neighborhood U of 7 go (r). Denote by r, r$ the r- 
components of integral curves 7, r ) qo of W. Then there are 5 > 0, a neighborhood U' 
of M", and a neighborhood V' of r so that for all < t\ < r and all integral curves 
7 of W lying in U' for the time interval (t%,t) and satisfying \r(ri) — ?*o( T i)| < 5 
we have that j(t') G U for all r' G V. This is because r' = F(r, tp, 9) where F is 
C 1 and the values of ip, 9 for 7 and 70 will be close together at all times in (ti,t) 
if [/' is chosen sufficiently small. 

Next, for the 5 > and neighborhood U' of M^, obtained above there is a 
neighborhood V of q m 2j so that for all q G V the curve 7 9 first runs inside 
{r < 6/2} and then inside U' . Now ^ = 1 + 0(r 2 ) implies that the travel time 
Ti(g) for the first part is at most on the order of 5/2. Then the r-components r q 
and ro of j g , j qo satisfy \r q (ri(q)) — ro(ri(g))| < <5, so the first part of the argument 
can be applied with 7 = 7 q and t\ — Ti(q). Summarizing, given any neighborhood 
U of 7o(t) we have found neighborhoods V of qo and V of r so that 7 g (r') G f/ 
for all q G V, r' G V' . This proves the continuity of exp aM , from the left with 
respect to q. Continuity from the right follows by the same argument applied to 
the left endpoint of Tj+i. The continuity of the inverse is proved in a similar way, 
but easier. For example, the second component, r, of exp^/ is simply the distance 
to the boundary, and its continuity follows from the triangle inequality. 



Remark 6.5. The proof also shows why the exponential map may be discontinuous 
in general: Suppose k > 2, i.e. the junction S has at least two maxima and at least 
two minima. Suppose T pi does not approach p' x but rather continues in the upper 
half plane 9 > and then approaches p' 2 . Examples of functions S yielding this 
boundary dynamics can easily be constructed. Let q\ , (72 be the labels of the geodesies 
leaving p[, p 2; respectively. Then, for any fixed r > 0, the point 7 g (r) will approach 
7g 2 (t) rather than j qi (r) as q — > q\ from the left, so exp^ M is discontinuous at 
(<ji,t) for any r > 0. Also, it is easy to see that this discontinuity cannot be 
removed by a simple reordering (i.e. by a different gluing prescription), which in 
any case would be unnatural since it would break up the order preserving property 
°fe x PdM- 
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Remark 6.6. The proof of Theorem \l.J\ gives more than is stated in the theorem. 
For example, the restriction of the smooth vector field V on T*M to the C 1 - 
submanifold M u may be projected to M . This yields a C l vector field V on M 
(which is C°° away from the finitely many geodesies leaving minima of S) whose 
integral curves are the geodesies starting atdM , as sets. The vector field W = r~ l V 
on M has, by definition, unit length everywhere, and it is the gradient field of the 
distance function to the boundary, W = Vdist. Note also that r' — £ + 0(r 2 ) 
1 + 0(r 2 ) yields that r — gdist for a non-vanishing function g equal to one at the 
boundary, and therefore we get V = gVdist 2 . So the gradient o/dist 2 has a certain 
amount of smoothness. 



Remark 6.7. The blowing-up procedure explained in Remark \ 6.2\ can be used to 
describe the boundary behavior of exp aM; in the setting of Theorem 1.4 Let M 



be the space obtained by blowing-up all m axi ma o f S in M with respect to the 



geodesic vector field V (see Remarks 6.2 and 6.6), with blowing-down map (5 



M — > M . Also, let M' be the space obtained from M by collapsing all intervals 
between adjacent front faces to points; these intervals correspond to intervals on 



dM from one maximum to the next. As explained in Remark 6.2 the S 1 part in the 
domain of expg M may be naturally understood as the union of all the front faces of 
(5, glued at their ends in cyclic order, hence with the boundary of M' . Then exp aM 
extends to a homeomorphism from S 1 x [0,tq) to a neighborhood of the boundary 
in M' . _ 

Note that M, M and M' are different ways of adding a 'boundary' to M. The 
correspondence of M' and M is analogous to birational maps in algebraic geometry. 
Both are obtained as blowings-down (contracting certain sides to points) of the 
manifold with corners M: 

M' < — M — > M 

Proof of Theorem |1.5[ Let p m j„ be a point where S has a local minimum and 
where S vip > |. Then the eigenvalues A2(p m in), A3(p m i n ) are non-real with real 
part — |, so nearby trajectories of V on T*dM spiral towards p min . 

Consider the maxima p_ , p+ of S closest to p m i n , where cp- < ip m in < l P+ for the 
corresponding parameters. W.l.o.g. we may assume S(tp-) < S(tp+). Then we are 



in the situation discussed before Lemma 5.3 (where we put </?_ = for simplicity), 
more precisely in the first case of (27 1. We claim that the trajectory T p considered 



there approaches p m j n as t — > oo. To show this, recall from Proposition 5.1 that the 
function S(<p) +8 2 /2 is strictly decreasing along T p . At t = —oo this function has 
the value S(ip—), and on the line tp = tp+ its values are at least S(ip+) > S(ip-). 
Hence, L p _ can not cross or approach this line, and neither the line ip = </?_. So 
the only possible limit point of L p _ for t — > oo is p m ; n , which was to be shown. 

Now note that the spiraling of L p _ around p min as t — > oo implies that the 
projection of L p _ to the t/3-axis oscillates around t/? m in infinitely often. 

Now consider the unstable manifold M p _ . The image of T p is an open subset of 
dM p _ . By continuity of V, for any T > and e > there are interior trajectories 
of V contained in M" that follow the boundary trajectory T p within an error 
of e on the time interval [0, T]. In particular, for any e > there is a trajectory 
7 starting at p_ whose projection to M intersects the projection of the trajectory 
Mp m . n at a point with r < e. This proves the theorem. 
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7. Examples 

We consider surfaces IcR 3 given as in ([2| with X = dX x R + a cylinder, for 
different boundary curves dX C M 2 . We use the Euclidean metric on M 3 . Write 
coordinates on R 2 as u = (v,w). Recall that S(v,w) = v 2 + w 2 , see Remark 



2.4 



(1) If dX is a circle centered at the origin then S is constant, the geodesies of 
X starting at the origin foliate the surface. This is obvious by rotational 
symmetry. 

(2) If dX is an ellipse centered at the origin, 

2 2 

— + — = 1, with c> b > 0. 
& b z 

The function S has two maxima at (±c, 0) and two minima at (0, ±6). The 
boundary dX is simply tangent to the circles centered at the origin passing 
through these points. Since dX bounds a strictly convex set containing the 
origin, Theorem |1.4| is applicable, see Proposition [53J There is one geodesic 
starting from each of the points (0, ±6), all the others start at (±c, 0). 

(3) We now consider a circle not having the origin in its interior, 

(v - c) 2 + w 2 = 1, c> 1. 

The function S has a minimum at p m ; n = (c — 1,0) and a maximum at 
Pmax = (c + 1,0). These are non-degenerate, so we have simple tangency 
again. In the arc length parametrization v(tp) = c + cos^, w((p) = sinip, 
we find S(ip) = c 2 + 1 + 2c cos if , so S vtp (ip) — — 2ccos ip. At the minimum 
ip = 7r we see that S w = 2c, which is also the maximal value of S vv . We 
obtain the following picture: 

There is a single geodesic 7 m i n starting at p m i n , all others start from 
Pmax- Their behavior depends on c: 



If 1 < c < 9/8 then Theorem 1.4 is applicable 



If c > 9/8 then we are in the situation of Theorem 1.5 Geodesies starting 
at p m ax almost tangentially to dX will approach p m m, then oscillate around 
7min many times before they escape far away from the singularity. 



Our first and simplest example above might suggest that Theorem |1.2| is not 
so interesting. However, the point of the theorem is that all that matters for the 
conclusion is the constancy of S on the boundary. Hence, any surface X with the 
same boundary will have a foliation by geodesies near the singularity. As another 
example, consider a surface X embedded in K" with n > 3, with dX any curve 
in a sphere R ■ S" -2 . Then there is no rotational symmetry, not even for dX, but 
Theorem |1.2| still gives a foliation by geodesies. 

8. Higher order cusps 

Everything carries over to higher order cusps, with minor numerical changes. 
Let k > 2. A cusp singularity of order k is one which can be resolved by an 
order k blowing-up (u, z) i— > (z k u,z), or alternatively by a sequence of k standard 
blowings-up, or alternatively by one standard blowing-up followed by an order 



(k — 1) blowing-up. Using this last characterization the proof of Proposition 2.3 



carries over, mutatis mutandis, and shows that the resolution of a cusp singularity 
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of order k yields a k-cusp metric on the manifold with boundary M — X, which is 
of the form 

g|l7 = (1 - k{k - l)r 2k - 2 S + 0{r 2k - r )) dr 2 + r 2k h 
where S and h are as before. The factor k(k — 1) is put in so that in the case of 



M = X the function S is determined by ( 11 ) as before 



The correct rescaling of the cotangent variable is 

Q = V 

and then the energy is given by ^£ 2 + r 2fe ~ 2 (Go + G) where 

Go = \k(k - l)S^ 2 + \C~ 1 9 2 . 
The rescaled geodesic vector field is V = rW and near the boundary yields the 



equations, generalizing ( 18 ) 



r = ri + ... £ = fc(l - £ 2 ) + • • • 

<p = (G )e + ... 9 = -(G Q ) V - (2k - l)£9 + 



where the dots indicate higher order terms in r. This has linear part as in ( 19 ) 
with S replaced by \k(k — 1)5 and 3Id by (2k — l)Id. Its eigenvalues (in the case 
of surfaces, say) are 1 and 

-(2k - 1) ± y/(2k - l) 2 - 2k(k - Yja 



where a = S vv in arc length parametrization. Therefore, the value of a which 
separates real from non-real eigenvalues is 

(2fc-l) 2 n 1 

«fe - TH-n = 2 -' 



2k(k-l) 2fe(A;-l) 

The analysis of the boundary dynamics carries over literally if one replaces S by 
i/c(fc — 1)5 and 30 by (2k — 1)9, and | by a&. The barrier function should be taken 

fc(fc-l) c 
(2fc-l) °f 

Since (2k — l)Ai(p) + A2(p) + A 3 (p) = 0, there are resonance relations among 
Ai (p), A2 (p) , A3 (p) . The value that the function S vv must avoid at a local minimum 
in the statement of Theorem |1.4| for fc-cusp is still 2, where ip is an arc-length 
parameterization of the boundary curve. Note that A2(p) = 1 — fc iff <z(p) = 2. The 



analogue of Proposition 4.6 holds, but when A2 ^ Q, we can locally linearize V by 
means of a G 2fc ^ 2 diffeomorphism. If A2 eQ\{l - fc}, we can only locally linearize 
V by means of a C 1 diffeomorphism. 

Now the proofs of all theorems carry over almost literally, with | replaced by 
afc. Note that strictly decreases to 2 as fc — ¥ 00, so the condition on S for having 
a foliation by geodesies near the singularity becomes (slightly) stronger the larger 
fc is. 

Finally, Theorem |1.4| applies to resolutions of cusp singularities of any order 
fc > 2 satisfying the convexity hypothesis of Proposition |5.5| since > 2 for all fc. 
This is precisely what is needed, see the statement of Proposition |5.5| 
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